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Abstract 



A linear spin wave analysis of dimerization of alternating Heisenberg system with spins s% and S2 on linear chain as 
well as square lattice is presented. Among the several possible dimerized configurations considered in two dimensions 
the plaqucttc configuration is found to be energetically the most favored one. Inclusion of a variable nearest neighbor 
■ exchange coupling J(a) = — leads to a uniform power law behavior: that is to say, the same (^-dependence is found 
(i) in chains as well as in square lattices; (ii) in systems consisting of different pairs of spins si and S2\ (iii) for the 
magnetic energy gain, the energy gap, the energy of the gapped magnetic excitation mode as well as for the sublattice 
magnetization; (iv) for all the configurations of the square lattice; and (v) in the entire range of 8 : (0 < 6 < 1). The 
variable exchange coupling also allows the energy of the gapped excitation spectrum to be 5-dependent even in the 
linear spin wave theory. 
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I. INTRODUCTION: 



Extensive interest is being currently shown in alternating spin systems consisting of two sublattices with unequal 
spin magnitudes s\ and S2 with a net non-zero spin per unit cell, as shown in Fig.l. Such systems are realized in 
bi-metallic chains with the general formula of ACu(p&aOH)(fi20)3.2H20 where pba OH is 2-hydroxy-l,3-propylenebis 
(oxamato) and A = Mn, Fe, C, Ni 0]. These ferrimagnetic chains are also referred to as alternating or mixed spin 
chains and are regarded as Heisenberg systems 



Alternating spin systems have been studied extensively by various techniques: by spin wave theory (SWT) [p|-^0[, 
spin wave expansion (SWE) |ll],[l^], Monte Carlo (MC) method pl| 10|, density matrix renormalization group 



(DMRG) technique |5 pC|Jl2| ], the method of matrix product (MP) states |3j] and by exact diagonalization (ED) 
method §|§. 

For an alternating-spin chain the zero temperature ground state energy and sublattice magnetization were evaluated 
using SWT SWE 0, DMRG [|§@ and QMC g. The results are summarized in Table 1. The linear 

spin wave theory gives higher values for the ground state energy and lower values for the sub-lattice magnetization 
compared to the more exact methods. Recently Ivanov et al. |ll],|l2| used a second-order spin wave expansion to 
calculate ground state energy as well as sublattice magnetization. Their results differ by 0.03% for the ground state 
energy and 0.2% for the sublattice magnetization from the DMRG results, as shown in Table 1. 

The thermal behavior was also investigated for ferrimagnetic chains |^|||] [l0) ■ Besides verifying the existence of two 
(gapped and gapless) excitation modes, the specific heat and magnetic susceptibility of ferrimagnetic chains were also 
shown to depend upon temperature as T 1 / 2 and T~ 2 respectively at low temperatures ||,[l0|]. It was also shown that 
this model behaved as a ferromagnet at low temperature, but as a gapped antiferromagnet at moderate temperatures. 

Modified spin wave theory, which includes Takahashi constraint, was also shown to give results in surprisingly good 
agreement with those from quantum Monte Carlo method in the thermodynamic limit of this system p|,|Ic|] . 

Dimerization of chains with spins s\ and S2 ( s% > s^) on alternating sites was recently |5|,^| studied the using the 
Hamiltonian 

H = [(1 + *) Sl,n ■ S 2 ,n + {l-5) S 2 ,„ • S hn+1 ] , (1) 

n 

where the total number of sites (or bonds) is 2N and the sum is over the total number of unit cells N. 5 is the 
dimerization parameter and is taken to vary between and 1. Linear spin wave theory and DMRG were used m 
to investigate the ground and low-lying excited states for both uniform and dimerized chains. In both the methods 
the ground state was found to be ferrimagnetic. One point of focus for us in the study of chains is that the LSW 
theory with the Hamiltonian in Eq.(Q) showed that the energy gap at k = in the gapped mode did not depend on 
the dimer parameter 5, while the DMRG predicted an almost linear dependence |^||.Th.e DMRG results on chains 
also show that the transition to a spin-Peierls state is conditional in that the ground state energy depends upon the 
dimerization parameter as 8 V with v=2 ±.01. 

This has motivated us to investigate a dimerized alternating spin Heisenberg model by using a linear spin wave 
theory using an ansatz of a variable nearest neighbor exchange coupling that was recently used to study dimerization 
in a uniform (single spin) Heisenberg system (r|]. We would also like to extend our work to square lattices. 

A need for describing nearest neighbor exchange interaction as dimerization sets in two dimensional lattices was 
recently discussed Jl3] | . Among the various lattice deformation modes which allow for dimerization, some require that 
the effect of the ensuing elongation or contraction of nearest neighbor distances be reflected in the nearest neighbor 
spin-spin exchange couplings. Since an exchange integral for a nearest neighbor distance a is roughly Jl4[ 

J(a) = -, (2) 
a 

we assume that when the nearest neighbor distances change from a to a(l ± S), the exchange couplings change from 
J to j^g. Thus, to linear order in 5, the interaction has the familiar form J(l ± 5). The form in Eq.(||) allows 
for incorporating changes in the nn exchange couplings in various situations of lattice deformations. It was shown 
||l3f that the logarithmic (^-dependence of various quantities like the gain in the magnetic ground state energy, etc., 
e g ~ rr-r i , both in one and two dimensional lattices, can also be a result of using this variable nearest neighbor spin-spin 
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exchange coupling. The form in Eq.(g) gives a logarithmic dependence over not just the S — > regime, but over the 
entire range of 5 from to 1. In what follows, we shall use for exchange interaction the form in Eq.(||). 

In this paper we will study alternating spin systems formed with different pairs of spin values; i, 1, and | using a 
zero temperature linear spin wave theory. We have considered three alternating spin systems from these spin values: 
(1, \) (denoting s\ = 1 and s 2 = \)\ (§, \) and (|, 1). We would like to see the effect of including the variable nearest 
neighbor exchange coupling on the (S-dependence of the physical quantities such as the gain in magnetic energy, the 
sub-lattice magnetization and energies of the excitation modes in both one and two space-dimensions and for different 
spin systems. We would also like to see if it gives a ((-dependence of the gapped excitation energy mode. In section two 
we will study these three alternating spin systems for a Heisenberg linear chain using the variable nearest neighbor 
exchange coupling. The energy and magnetization of such systems will be computed using LSW theory. Critical 
exponents of the dimer alternating chains will also be calculated. We shall then study alternating spin systems on a 
square lattice for several proposed dimer configurations in section III. 

The question of frustration in a ferrimagnetic chain or a ladder due to an antiferromagnetic second neighbor 
interaction has also received some attention recently. It was shown that a strong frustration leads to disappearance of 
the long range ferrimagnetic order through a discontinuous transition to a singlet state [^2| , p^|Jl^ 1 . It was also shown 
that the spin wave theory can predict realistic results for a frustrated system at least for the case of weak frustration 
|l2|| . We shall attempt to see the effect of the variable nearest neighbor exchange interaction on weakly frustrated 
chains and square lattices in a future publication. 



II. ONE DIMENSIONAL ALTERNATING SYSTEM 



The alternating dimer Hamiltonian on a chain with two spins s\ and s 2 can be rewritten using the variable nearest 
neighbor exchange coupling defined in Eq.(S) as 



H = J2 



Sl,2i ' $2,21+1 + ~ ^2,2';+! ' "5*1,24+2 



1 



(3) 



A linear spin- wave analysis is usually performed with the help of Holstein-Primakoff (HP) transformations to bosonic 
spin-deviation operators. For the two sublattices the HP transformations are: for spin-si 



°l.n 


= (2si 




= 4(2 


°l,n 


= Si - 



and for the second sublattice with spin-S2 



9+ - 



5'. 



2,n 



si - a^a„) 1/2 



6t(2 S2 -6t6„) 1 / 2 
(2s 2 - b\b n ) x l 2 b n 



S 2,n = b l b n ~ S2 



(4a) 
(4b) 
(4c) 

(5a) 
(5b) 
(5c) 



where Si is the magnitude of the spin on sublattice i. A linearized Hamiltonian is obtained by substituting HP 
transformations into Eq.(|^) and keeping terms up to the quadratic order in the spin-deviation operators a and b. The 
linearized Hamiltonian in Fourier transformed variables is 



H = [Aia\a k + A 2 b\b k + B(k) (4&j. + Mfc) + C 

k 



(0) 



with 



Ai = J p s 2 
A 2 = JpSi 
B(k) = A kA /I^2~ 

C = —J P S\S2- 



(7a) 
(7b) 
(7c) 
(7d) 
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Here 



and 



A fe = ^(J p cos(fc)) 2 + (J m sin(fc)) 2 (8) 



Jrn = JpS (10) 



to 



The linearized Hamiltonian in Eq.(^) can be diagonalized using Bogoliubov transformations 

a k = u k a k + v k 0[ (11a) 
b k = u k (3 k + v k a\ (lib) 

h = Y,[ £ 9 + e i + Mwlfik] (i 2 ) 

k 

where the coefficients u k , v k are constrained by the condition ui — vt = 1, a k and /3 k are the normal mode boson 
operators, E\(k) and E 2 (k) are the energies of the two excitation modes and e g is the ground state energy per site, 
with u(k) and v(k) defined as 



< k ) = \ \7 (13a) 





f A 2 


+ 6= 






\-A 2 


-6 



"*> = V £ (13b> 



& = \/(/li + il2) 2 -4B(t) 2 (13c) 

The two excitation modes are 

E 1 (fc) = (^ 1 -A 2 +6)/2 (14a) 

£ 2 (fc) = (^ 2 -A 1 +6)/2 (14b) 

It is easy to see that E\ is the gapless mode and E 2 has a gap. 
The ground state energy per site e g is given by 

s g = C-A 1 -A 2 + J2& (15) 

k 

and the staggered magnetization in the two sub-lattices corresponding to the spins s% and s 2 respectively is 

Mi = Si- < L> > (16a) 
M 2 =< D > -S 2 (16b) 

with < D >=< a\cii >~< bjbj > is the average taken in the ground state, which is the Neel state, at zero temperature. 
This average of spin deviation operators can be determine from 



N 

k 



with k running over half the Brillouin zone. 
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For the three kinds of alternating spin chains, referred to as (1,4), (§, §) and (|, 1), the ground state energy, the 
excitation energies and magnetization can now be calculated as functions of the dimerization parameter 8 . Previous 
calculations invariably took spin-spin exchange couplings alternately as J(l ± 8), which, as mentioned above, can be 
taken as an expansion of the interaction in Eq.(^J) to the first order in 8, implying that the results are valid only in 
the critical regime <5 — > 0. The advantage of taking the variable nearest neighbor exchange coupling is that the results 
will then be valid also in the limit 8 — > 1. 

The ground state energies e g per site for undimerized chains, 8 = 0, were found, expectedly, to be the same as those 
found earlier H for the three systems. After including dimerization, our calculations also confirm that the ground 
state energy of all the three systems described above decreases with 8. This is shown in Figure 2, where energy gain 
s g (8) — e g (0) is plotted against 8. Numerical fitting shows that, as against earlier results, the magnetic energy gain 
Eg (8) — £ g (0) has a logarithmic dependence on 8, t^jt, for the three systems discussed here, with values of v between 
1.4— 1.6 in the entire range < 8 < 1. Fig. 2 shows that the chain (|, 1) has higher gain than the other two systems. 

As expected, our calculations also find two branches of the excitation spectrum, one gapless and the other with a 
gap at k — 0, in the three systems. 

As stated above, the LSW theory with the spin-spin exchange coupling J(l ± 8) in the presence of dimerization 
allowed no 8— dependence of the energy gap in the second mode, while the DMRG found almost a linear 8— dependence 
for these spin systems ]^,^| . We find that by including a variable nearest neighbor coupling constant defined in Eq. 
the spin wave theory also allows for a 8— dependent energy gap A(<5) = E2(8) — £ g (8) in the second mode. This is 
because the <5-dependent terms do not now cancel out for the gapped excitation mode at k = as they did with the 
coupling J(l± 8). The dependence is found to follow the same logarithmic behavior, r^-jr, as the ground state energy, 
with v varying between 1.4 — 1.6 for the entire range < 8 < 1. This is true for all the three spin systems defined 
here, and is larger for the (|, 1) system than the other two (see figure 3). 

The staggered magnetization M{8) was also found to follow the logarithmic (5-dependence up to 8 < 0.5, but follows 
a different behavior for 8 > 0.5. Mi (8) against dimerization is shown in Fig. (4) for the three systems. Again we see 
from Fig. (4) that the chain with (|, 1) has higher value of magnetization than the other two systems. 

It is worth mentioning here that by using coupled cluster method p3|, we had found that the ground state en- 
ergy and the staggered magnetization of a spin-half Heisenberg chain follow the same logarithmic behavior using the 
variable exchange coupling defined in Eq.(|^), in both small and large values of dimerization 8. This gives us more 
confidence about the results we have obtained by LSW theory. 



III. TWO DIMENSIONAL ALTERNATING HAMILTONIAN :- 



The dimerization on two dimensions lattices differs from that on chains. There are several ways in which distortions 
of a square lattice can occur, each one of the possible configurations giving a different dependence of the ground state 
energy on the dimerization parameter fl3|| . 

We will use some of these configurations, illustrated in Fig. (5), to study the alternating spin square lattices. The 
alternating dimerized Hamiltonian for a two dimensional system can be written in general as 

H = ^2 [Jx,fj,Sl,i,j ' S2,i+fi,j + Jy,fiSl,i,j ' ^2,iJ+fj] (18) 

where the indices 1 and 2 on the spin vectors refer to the two sublattices with spins of magnitude s± and 82- For reasons 
described earlier, we use variable nearest neighbor exchange couplings. These are defined for different configurations 
as follows: 

Configuration (a) 
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Configuration (b) 



Configuration (c) 

Jx.fi = Jy,p = (1+fj.S) — ~ t 1 ^ 

Configuration (d) 

Jy u — ~, i , £Si J ( 1 - fiS ~ (1 - 4)t 2 ) 

Configuration (e) 

J^ = J v ^ 7W ^^j(l-,8-il-^) 

We would like to investigate the five configurations resulting from dimerization of a square lattice in order to see (i) 
which one of these leads to the largest gain in magnetic energy as the dimerization sets in, (ii) if the use of variable 
exchange coupling leads to a single power law behavior valid for the entire range of 5, and how the law differs from 
that in the case of chains, (iii) the (5-dependence of the second mode of excitation E2, (iv) the behavior of staggered 
magnetization, and (v) the generality of these investigations regarding the three spin systems discussed here. 



The linear spin wave analysis follows the same procedure as for the chain above. The same equations are applicable 
in this case, but the various coefficients entering the theory have now the following values: 

Ax = J p s 2 (19a) 

A 2 = Jpsx (19b) 

B(k) = r(fc)VsiS2 (19c) 

C = -J p sxs 2 (19d) 

where 



T(k) = \J (J px cos(k x ) + J py cos(ky)) 2 + (J mx sin(k x ) + J my sin(k y )) 
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and 

Jp = (Jx,+1 + Jx,-l + Jy,+1 + Jy,-l)/^i 
Jpx — (Jx, + 1 + Jx,-l)/4, 
Jpy = {Jy,+1 + «/y,-l)/4, 
Jmx — (Jx, + 1 Jx, — 1)/4, 
Jmy = {Jy, + X — Jy, — 1)/4- 



The ground state energy e g (5) defined in Eq.(h5h, energies of the two excitation modes Ei(k) in Eq.(14) and stag- 



gered magnetization M($) defined in Eqs.([19|) can now be calculated as functions of the dimerization parameter 6. 

The ground state energy e g (S = 0) is found to be —1.2,-1.7158 and —3.3709 for the three spin systems (1, |), 
(|, i) and (|, 1) respectively. Staggered magnetization Mx{M2\ on the first {second} sublattice is 0.8907 {—0.3907}, 
1.4241 {-0.4241} and 1.3597 {-0.8597} for the three systems . These values are listed in Table 2. 



Our calculations confirm that, like in chain, the gain in magnetic energy increases with 5 in all the proposed config- 
urations. This is shown in Figure 6, where the energy gain s g (S) — e g (0) is plotted against 6 for the five configurations. 
It also shows that the plaquette configuration of Fig. 5(c) is energetically the most favorable state, while there is 
hardly a discernible difference among the configurations (a), (b) and (d). It is also interesting to note that the mag- 
netic energy gain under dimerization of an alternating spin square lattice also varies as r^jr with u = 1.4— 1.6 in the 
entire range 0< 6 < 1, exactly as in the case of a chain. This is singularly an effect of taking the variable exchange 
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coupling denned above. 

The S dependence of the energy gap, A(<5) = E-2(5) — e g (S), for the five configurations is shown in Fig. 7, showing 
greater stabilization of the dimerized state with increasing 8. We also find that, like the magnetic energy gain, the 
energy gap increases with 8 as J^gr in the small 8 regime for all the five configurations with v = 1.4 — 1.6. The 

configurations (a) - (d) also have the same dependence on 8 in the entire range of 8. The difference between the 
dimerization of a square lattice for these configurations is again markedly brought out in Fig. 7. 

Our calculations give staggered magnetization for the un-dimcrized alternating spin square lattice M\{8 = 0) = 
0.8907, 1.4241 and 1.3597 and M 2 (8 = 0) = -0.3907, -0.4241 and -0.8597 for (1,±), (§,±) and (|,1) respectively. 
As dimerization sets in, magnetization decreases in all the configurations we have chosen, as shown in Fig. 8. This is 
also the case for the entire range < 8 < 1, except in the case of configuration (e) for which the magnetization rises 
again after 8 > ^. 

Configuration (e) is peculiar in the sense that 5 — ^ is a special point for it; the shorter bond length is symmetric 
about this point, having a minimum value of -j=. At this point the distortions give rise to a rectangular lattice with 

sides \/2 and The energy gain increases with 8 up to S = j, and then decreases. 

For all the five configurations, we found that the magnetization also varies as t^j-\ in the small 8 regime with the 
exponent v = 1.4— 1.6, exactly as the energy gain and the energy gap. However, while for configurations (a-d) in the 
full range < 8 < 1 the magnetization follows the same power law with the exponents as v = 1.4 — 1.6, configuration 
(e) has a distinctly different behavior in this regime. 

In summary, we have studied the spin-Peierls dimerization of an alternating spin Hcisenbcrg system on a chain 
and a square lattice taking variable exchange couplings based on the ansatz J(a) = ^, for three kind of alternating 
spin systems, namely (1, \), (§, \) and (|, 1). We have included different possibilities of dimerization in the case of 
a square lattice. The ground state energy as well as staggered magnetization decrease continuously with increasing 
dimerization in both ID and 2D . In 2D, the plaquette configuration with dimerization taking place simultaneously 
along both the principal square axes has markedly lower ground state energy and magnetization than the other config- 
urations; (a), (b), (d) and (e). The plaquette configuration stands out as the most favored mode of dimerization. The 
energy gap also corroborates the above conclusions. It has also been shown that the magnetic energy gain, energy gap 
and staggered magnetization follow a uniform dependence upon the dimerization parameter 8 as (i) in chains as 

well as in square lattices; (ii) in systems consisting of different pairs of spins s\ and S2] (iii) for the magnetic energy 
gain, the energy gap, the energy of the gapped magnetic excitation mode as well as for the sublattice magnetization; 
(iv) for all the configurations of the square lattice; and (v) in the entire range < 5 < 1. The variable exchange 
coupling also allows the energy of the gapped excitation spectrum to be <5-dependent even in the linear spin wave 
theory. 

We would like to thank to Dr. D. Sen for providing Refs. 5 and 6. 
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Figure captions 

Figure 1: A schematic sketch of an alternating spin chain. The larger and smaller arrows indicate the larger (si) 
and the smaller (s 2 ) spins. The hollow (filled) circles represent the positions of spins in the undisturbed (dimerized) 
chain. 

Figure 2: The magnetic energy gain s g (S) — s g (0) vs the dimerization parameter 5 for ID alternating spin chain in 
the full range of the dimerization parameter < 5 < 1. 

Figure 3: (^-dependence of the energy gap as dimerization sets in for different alternating spin chains. 

Figure 4: (5-dependence of the staggered magnetization for one of the two sublattices, Si, for the three alternating 
spin chains. 

Figure 5: Five configurations for the dimerization of a square lattice, (a) a columnar configuration caused by a 
longitudinal (tt, 0) static phonon; (b) a staggered configuration caused by a (tt, tt) static phonon with polarization 
along .x-dircction. Like (a) , the dimerization occurs along one direction only, but the sequence of alternate couplings 
itself alternates along the other direction; (c) dimerization along both the directions, caused by (tt, 0) and (0,7r) 
phonons, making a plaquette of four nearest neighbour spins; (d) again, dimerization along both the directions, but 
taken staggered along the vertical direction; (e) another staggered dimerization that is caused by a longitudinal (tt, tt) 
phonon. Large arrow belongs to the first sublattice and short one belongs to second sublattice. And the open circles 
indicate the square lattice sites and the solid ones show the dimerized lattice. 

Figure 6: The gain in magnetic energy e g (S) — s g (0) as dimerization sets in with increasing S for the five configu- 
rations of a square lattice in the range < 5 < 1 for (a) spin (1, 5), (b) spin (5,5) and (c) spin (|, 1). 

Figure 7: Dependence of the energy gap A on 8 for the five dimerization configurations of the alternating square 
lattices for (a) spin (1, 5), (b) spin (5,5) and (c) spin (|, 1). 

Figure 8: (5-dependence of the staggered magnetization of an alternating spin square lattice calculated for the five 
dimerization configurations, (a) spin (f , 5 ), (b) spin (|, 5) and (c) spin (|, 1). 
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Table 1: Summary of the ground state energy per site and sublattice magnetization values calculated by different 
methods for alternating spin chains made up of the three spin systems namely (1, 5), (§15) ancl (§> !)• 



Spin 


Method 


c 9 


Mi 


M> 




Mr _ 


— 0/245 


r\ 770 
0. (79 


— U.Z79 




— u. < / /o 


O 70Q 

U. ( »o 


— U.zyo 


T CJWT I 1 
LO VV 1 4_ J 


j l n 71 r 


u.oy 


1 Q 1 ^ 


JJIVIxILl . 




U. i y4zo 


n ooO/i s 
— U.zyz4o 


SWE 4 1 


1 -0 72715 


0.79388 




MSW kl 


-0.7295 






(» I) 


LSWT i] 


-0.979 


1.315 


-0.314 


DMRGL^ 


L -0.98362 


1.35742 


-0.35742 


SWE iiS'L. 


1 -0.9834 


1.3666 




(1-1) 


LSWT 1 3 


-1.914 


1.040 


-0.540 


DMRG.^ 


I -1.93096 


1.14427 


-0.644 


SWE 


1 -1.9316 


1.1461 





Table 2: The ground state energy per site and the staggered magnetization of the un-dimerized alternating spin 
square lattice for the three spin systems as calculated in the linear spin wave theory. 



Spin system 




Mi 


M 2 


(U) 


-1.2 


0.8907 


-0.3907 


r± I) 


-1.7158 


1.4241 


-0.4241 




-3.3709 


1.3597 


-0.83597 
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Fig. 6 (a) 
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Fig. 6(b) 
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Fig. 6(c) 
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For Spin(3/2,1) 
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Fig. 7(a) 
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Fig. 7(b) 
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Fig. 7(c) 
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Fig. 8(a) 
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